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The broad application range of the predator-prey modelling enabled us to apply it to 
represent the dynamics of the work-employment system. For the adopted period, we 
conclude that this dynamics is chaotic in the beginning of the time series and tends to 
less perturbed states, as time goes by, due to public policies and hidden intrinsic system 
features. Basic Lotka-Voltcrra approach was revised and adapted to the reality of the 
study. The final aim is to provide managers with generalized theoretical elements that 
allow to a more accurate understanding of the behavior of the work-employment system. 
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"Economies possess general ecosystem properties, such as dynamism, evolution, integrity, 
stability and resilience. Economies are inextricably embedded in larger natural ecosystems, and 

exchange flows of materials and energy with natural systems" . 



S. Farber & D. Bradley 



1. Introduction 

The concept of equilibrium in the predator-prey population dynamics has its ori- 
gins in the works of Lotka ^ and Volterra making the base of several theoretical 
models of interaction among species, with applications in the con text of ecological 
systems and wildlife management (Caughley and Sinclair, 1994^21 Earn et al., 2000 
231 Blasius et al, 1999^, including fluctuating environments (Collie and Spencer, 
19970% Particularly, Blasius et al. were very benefited by previous works on phase 



1 



January 12, 2013 22:45 WSPC/INSTRUCTION FILE Econophysics 



2 Nilo Sylvio Costa Serpa 

synchronization phenomena in coupled chaotic systems (Rosenblum et al., 1996 ^j). 
Also, there are generalized approaches of the Lotka-Volterra model (Tu and Wilman, 
1992^3). There is special interest on the problem of the growth of two populations 
conflicting with one another, known as "the problem of Volterra" The generality 
of type predator-prey models becomes possible to abstract them from their early 
ecological roots, and, by analogy, to apply them on a large range of mathematical 
modelling problems. Recently it befell to us to investigate the work-employment 
dynamics in Brazil by the implementation of a type predator-prey model. It is not 
the case of a classical application of the standard Lotka-Volterra competition mod- 
elling, since the variables storage dissimilarities among population registers, not the 
populations themselves; this is so because the existing cyclic pattern which entan- 
gles the work-employment evolution is much more easily discerned by dissimilar 
modelling. Besides, while is tempting to discard the standard Lotka-Volterra model 
as too simplistic, there is a real situation where the detailed and complex dissimilar 
modelling holds undeniable utility. The aim of this article is to show that the work- 
employment system in Brazil admits a predator-prey modelling. The model started 
from a simple correlation suggested by the superposition of two time series, one to 
the number of employed workers, and other to the number of active employers. The 
entanglement of the two series is analyzed with the aid of the clustering method- 
ology. In the present article, the abstract analogy with a natural ecosystem is the 
work market; the employer or corporation plays the rule of the predator; the prey 
is represented by the worker. The assumption of employers as predators has many 
reasons: 

• The number of employers is much smaller than the number of workers (as 
in linx-hare, lion-zebra, crocodile-gnu and cheetah-baboon predator-prey 
relationships). 

• The employers handle the number of workers depending on the profit mar- 
gins which they want, heavily controlling the levels of unemployment. 

• Employers often fail to collect the so-called "Guarantee Fund for Time 
in Service" related to their workers, causing irregularities that hinder the 
realization of labor rights in the event of contract termination agreement. 
Moreover, even without a contract being rescinded, the nonpayment of the 
Guarantee Fund blocks the employee to request it to buy a home. 

• Employers often fail to inform the CAGED (Cadastro Geral de Empregados 
e Desempregados jf| on the movement of their staff, resulting in bureaucratic 
disarray that blocks the access to some labor rights in the event of loss of 
working papers; furthermore, the fail to inform the CAGED affects the 
labor statistics. 

• Many employers absorb the workforce without formal contracts, depriving 
workers of their entitlements. 



a CAGED is a great database updated monthly and containing all workers and employers in Brazil. 
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• However the scenario is changing in Brazil, many employers still do not 
make investments in quality, training employees and providing opportuni- 
ties for better wages by productivity and capacity. The result is a high turn- 
over, low wages and a working class that never achieves a status of partic- 
ipation on the employer's profits. In sum, we have a highly non-egalitarian 
society which does not prioritize education and manpower quality. Socially 
speaking, the predator is who foment social discrepancies and, in present 
case, they are employers (public or private). That is the point here. 

From this point of view, we have not a model based on political preferences, as some 
watchers would think, but based on hard facts. There is no influence of maniqueist 
ideas, but simply findings. Paraphrasing Lawson, the offer of labor posts depends 
on the investment decisions of capitalists (employers) , a nd t he investment decisions 
of capitalists depend on the existing offer of manpower ^^J. Thus, since employers 
impose financial directions on the work market, often indifferent to the zeal for 
the rights of the working class, thus creating distortions and fueling widespread 
dissatisfaction, it is reasonable to consider the adopted viewpoint in constructing 
the model. Obviously, in a more egalitarian society, firmly established on ethics and 
education, economically viable and with eyes in the future, it would be required 
a comprehensive review of the basic premises assumed. Also, in face of the facts 
itemized above, employers appear in many situations to dampen fluctuations in 
employee populations. 



2. Observational premises 

The data sources were the CAGED, and the RAIS (Relagao Anual de Informagoes 
Sociais), two of the more important corporative databases of the Brazilian Min- 
istry of Work and Employment. First of all we looked for harmonic patterns in the 
behavior of the working class and active employers. Figure 1 clearly presents the 
referred pattern in the chosen period (1996-2008). We defined active employer or 
company as one that presents variation in the number of admitted workers from 
a certain month to another. Secondly, we verify the presence of lags between the 
harmonic patterns of both series. This feature is also clearly visible in figure 1, with 
relaxation of the system all along the last fifty months. Table 1 shows a segment of 
the CAGED for the first year of the time series considered in present study 
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Figure 1: harmonic patterns from workers (full line) and active employers (dashed line). 
The phase difference is clear in most of the peaks and valleys, especially in the fifty more 
recent months. 



Table 1: aggregate data 



Period 


Balance workers 


01/1996 


-12.626 


02/1996 


-4.094 


03/1996 


10.003 


04/1996 


118.918 


05/1996 


172.930 


06/1996 


115.028 


07/1996 


68.920 


08/1996 


46.937 


09/1996 


88.964 


10/1996 


19.466 


11/1996 


-15.899 


12/1996 


-258.516 



of the CAGED database. 



Workers Active employers 



23.743.110 


336.946 


23.739.016 


334.349 


23.749.019 


359.248 


23.867.937 


349.273 


24.040.867 


344.896 


24.155.895 


326.832 


24.224.815 


338.098 


24.271.752 


339.688 


24.360.716 


327.525 


24.380.182 


329.672 


24.364.283 


308.009 


24.105.767 


298.644 
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Figure 1-a: the relationship of workers, active employers and balance of workers. 



The original Lotka-Volterra predator-prey model makes several simplifying assump- 
tions; it was adapted in many ways, as we may see in the vast literature on predator- 
prey interactions from which the authors selected some of the more important pub- 
lications (see references). In many cases the former assumptions are relaxed or 
rearranged to fit some particular dynamics. There is no doubt that the predator- 
prey model adapts, at least conceptually, the work-employment scenario, as it is 
very reasonable to admit that employers and workers are, after all, two populations 
conflicting with one another. Since the formulations of Karl Marx - with the anal- 
ysis of the conflicts between worker and employer, and of the entailments of such 
conflicts with structural elements formed by the connections among politic, social 
and economic plains - we understand that conflicts belong to the dialectic nature of 
the capitalist work relations. Also Richard Goodwin, in a different approach, found 
in Marx arguments conceptual similarities that lead to his predator-prey dynamic 
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3. Methods and simulations 
3.1. Formal representation 

The pertinent information such a model can render is whether or not population 
abundances tend to an equilibrium at which both, workers and employers, will coex- 
ist and survive with few conflicts. The adopted model was based on the approach of 
Blasius and Stone ^ about oscillations with Uniform Phase Evolution and Chaotic 
Amplitudes (UPC A). An ecological UPC A model was chosen because it is suitable 
for treatment of two correlated periodic phenomena with chaotic amplitude varying 
in time, however, almost constant amplitude-independent frequencies as hinted in 
Figure 1. We applied the model to investigate the complex temporal phase synchro- 
nization in work-employment system at Brazil. In other words, our approach was 
based on the theory of synchronization of chaotic oscillations, defined in the most 
general case ^ as the entanglement between the phases of two coupled systems, 
while the amplitudes stay over chaotic regime in time. The real entanglement of 
phases in predator-prey relationships is much more complex than one could sup- 
pose at first. As pointed out by Holling, predator-prey interactions have shown that, 
in some cases of population densities, predators may control the number of preys, 
but this is not true to all density cases SI. However in Figure 1 populations grow 
and fluctuate in a net periodic way, the tests on the tangle between these popula- 
tions goes far beyond the simplicity of that plot, as the reader can see in Figure 
1-a. The balance of workers is a finite difference that indicates whether the work- 
employment scenario is favorable or not to the employee, and to what degree; it 
reflects the potential of the Mean (work market) to sustain the system evolving, by 
analogy with the available natural resources in an ecological niche. In other words, 
if the balance is negative we believe that there is a retraction of work market to sup- 
port occupied posts and, virtually, active employers. It is a variable that provides 
additional qualitative (positive or negative sign) information. Figure 1-a shows the 
active employers-workers-balance of workers relationship, illustrating the potential 
oscillations in these populations. The Lotka-Volterra-Serpa differential equations 
consider a state vector formed by the product between a geometric tracing vector 
and a dissimilarity vector of observables. We use pairwise dissimilarities between 
registers instead of the quantities directly stored in data base just to compute the 
relative differences which are really cyclic all along the time series. That is the basic 
difference from the classical Lotka-Voltcrra approach and the name was given by 
fellow researchers in honor of the author, Nilo Sylvio Costa Serpa. The dissimilarity 
matrix was made from the monthly populations of active employers and workers, 
and from the monthly balance of workers. This matrix computes all the pairwise dis- 
similarities (or " distances") between observations in the CAGED data set by means 
of the DAISY algorithm that executes mixed measurements. We assume that 
the pairwise differences computed from the time registers determine the dynamical 
configuration of the system. Thereby, our model has three differential equations, one 
for the dissimilarity of balance of workers u (resources) , one for the dissimilarity of 
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Figure 2: UPCA oscillations of employer (grey) and worker (black) dissimilarities with ini- 
tial time-delayed growth (a 2 = 1.0). 

workers v (preys) and one for the dissimilarity of active employers w (predators), 
such as 

' u = a.u.u — a.\.f\(u.u , v.v a , k\) 

v= —b.v.v o + 0ii.fi(u.u o ,v.v o ,ki)— 

-a 2 .f 2 (v.v ,w.w ,k 2 ) 
k w= - c.(w.w - w^) + ot2-h{v.v , w.w , fc 2 ), 

where /i and f 2 represent either the Lotka-Volterra term fi(x,y) — xy or the 
Holling type II term fi(x,y) = xy/(l + hx), and denotes the minimum level of 
existent employers dissimilarity when there is scarcity of workpower or doldrums. 
The overdot indicates differentiation with respect to time. The quantities a, b and 
c are growth rates of balance of workers, workers, and employers respectively. The 
Holling type II functional response SI introduces a decelerating intake rate related to 
the assumption that the dissimilarities of active employers is limited by its capacity 



(1) 
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Figure 3: UPCA oscillations of employer (grey) and worker (black) dissimilarities with ini- 
tial time-delayed growth (a 2 = 1A). 

to regulate dissimilarities of workers. Comparing to the original equations applied 
by Blasius, 

u = a.u — v, k\) 

v = - b.v + cti.fi(u,v, fci) - a 2 -f2(v,w,k 2 ) 
w= - c.(w - u> t ) + a&./aCv, w, k 2 ), 

the reader can see that we merge variables u, v, z with uq,Vo,zq and this changes 
the final geometric form of the system integration according to de observed data. 
Figure 5-b shows an example of UPCA oscillations for the last system. In our 
model, the positive constant k 2 measure the carrying capacity of employers to absorb 
manpower; the positive constant ki measure the carrying capacity of the work 
market to accept new workers. As pointed out by Kenneth Arrow and colleagues^, 
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Figure 4: UPCA oscillations of employer (grey) and worker (black) dissimilarities with ini- 
tial time-delayed growth (cti = 1.8/ 



" Carrying capacities in nature are not fixed, static, or simple relations. They are 
contingent on technology, preferences, and the structure of production and consump- 
tion. They are also contingent on the everchanging state of interactions between 
the physical and biotic environment. A single number for human carrying capac- 
ity would be meaningless because the consequences of both human innovation and 
biological evolution are inherently unknowable 1 ' . 

So, it is interesting to investigate what would happen for work-employment scenar- 
ios with very low carrying capacity of employers to create new posts. For instance, if 
we take fc 2 = we means that the employers are saturated or temporarily disabled 
to open new posts, although k\ remains greater than zero, that is, the work market 
remains in flux, but mainly due to high turn over. The contents of the differential 
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Figure 5: UPCA oscillations of employer (grey) and worker (black) dissimilarities with ini- 
tial time-delayed growth (oti = 2.0). 



variables were defined by means of the vectors, 

V (u < -w[l],v < -u[2],w < -u[3}) , 

U (u < -Tl[S, l],v < -U[S, 2},w < -U[5, 3]) , [ ' 

being ( j in [1, 2, 3]) the geometric components describing the state trajectories 
and U[S, 2j] the j dissimilarity components between observations with S registers 
calculated according to the definition of Euclidian distance, that is, 

( d r\ ^ 

n[5,i] oc I ^2 \ x Sk - x jK \ r j , (3) 



where d is the phase space dimensionality; x$ K and xj K are, respectively, the K — th 
components of the S — th and j — th registers. 
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Figure 5-a: UPC A oscillations of employer (grey) and worker (black) dissimilarities with 
initial time-delayed growth and &2 = 0.00008 (ct2 = 1.8/ 



3.2. The equilibrium condition 

No changes with respect to time in populations and resources, or in their dissimi- 
larities, are represented by setting the equations equal to zero. So, the equilibrium 
condition for the system (1) requires, 

a.u.u — 0ii.f\(u.u o , v.v , fci) = 
-b.v.v + ai.fi(u.u Q , v.v , fci)- 
-a 2 .f2(v.v 0l w.w 0l k 2 ) = 
-c.(w.w - + a2-f2(v.v ,w.w ,k 2 ) = 

Introducing the term of Lotka-Volterra (the simplest one) , the equilibrium - where 
the dissimilarities of active employers, resources and workers do not vary - leads to, 

a.u.u — cti.u .v .u.v = 

—b.v.Vo + ot\.u .v .u.v — a2.v .w .v.w = (5) 
— c.(w.w — w^) + a2.v .w .v.w = 



(4) 
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From the first equation we obtain, whereas population dissimilarity shall be never 
equal to zero, 

u.(a.u — ai.u .v .v) = 0, (6) 
a.u a — a\.u .v .v = Q, (7) 

(8) 



ai.u .v ai.Vo 
Taking the third equation we get, 

— C.W.Wo + C.W^ + Ot2-V .W .V.W = 0, (9) 

w (~c.w + a2-v .w .v) — ~c.w\ (10) 

(11) 



—c.w^ c.w^ 



- (c.w - a 2 .v .w .v) c.w - a 2 -v .w .v 
Finally, taking the second equation it follows, 

ai.u .v .u.v=a2-v .w .v.w + b.v.v , (12) 

a 2 -v .w .v.w + b.v.v 

u= , (13) 

a\.u .v .v 



a2-w -w + b 

u = . 

a\.u 

Substituting the value of w given by equality (11), it comes, 



(14) 



i c - w ,| 1 
u=\a 2 -w . ho , (15) 

C.W - C(2-V .W .V I a\.u 
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u=\ + b\ , (16) 

\c — a 2 .v .v J a\.u 

_ a 2 .c.w^ +b.(c- a 2 .v .v) 

U — - r . (17) 

a\.u . (c — a 2 .v .v) 

Now, applying the result set by expression (8), we have, 

t . , ( ai.c-a 2 .a\ 1 

u = a 2 .c.w< +b. [ -, (18) 

\ ai J u a . (ai.c - a 2 .a) 

ax-u 2 .c.w^ + b. (ai.c - a 2 .a) 

U= r . (19) 

a\.u . (ai.c — a 2 .a) 



For the term of Holling type II, that is, 

/i = u.u .v.v / (1 + ki.u.uo) , 

f 2 = V.Vq.W.Wq/ (1 + k 2 .v.v ) , 

the equilibrium sets, 

' a.u.u — ol\.u.uq.v.Vo/ (1 + k\.u.Uo) = 
-b.v.v + ai.u.UQ.v.vo/ (1 + ki.u.uo) — 

-a 2 . V.Vq.W.Wq/ (1 + k 2 .V.V ) = 

-c.(w.w — w^) + a 2 .v.v .w.wo/ (1 +k 2 .v.v ) = 

The system becomes much more complex. Nevertheless, we fix k 2 = 
that the employers are temporarily disabled to open new posts, which 
the solution as it follows. The second equation leads to, 

a\.u.uo/ (1 + ki.u.uo) — b 



w = 



a 2 .wo 



The first equation gives, 



a = ai.v.vo/ (1 + k\.u.uo) , 
a. (1 + ki.u.uo) = ai.v.vo, 
1 + ki.u.u a = 
ai.v.vo 



a\.v.vo 



a 
a 



UQ.a.ki 

Finally, the third equation sets v as, 

a 2 .v.v .w.w / (l +k 2 .v.v a ) = c.(w.w — w^), 

c.(w.w — w') 



v = 



o> 2 .vq.w.wo 



(20) 



(21) 



0, meaning 
determines 



(22) 

(23) 
(24) 
(25) 

(26) 

(27) 
(28) 
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Figure 6: the system in three-dimensional phase space. 




Figure 7: another vision of the same system. 

4. Discussion 

Insecurity generated by the employers and the work market, over many years of 
restructuring and adjustment to globalization as the Nineties, led to an exacerba- 
tion of individualism. Competition became more dominant than cooperation. Fear 
and instability certainly promoted higher levels of depression, while millions of peo- 
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Figure 8: workers and employers in phase space with an attracting central curve. 

pie reached the productive age, pressing more and more the formal work market. 
Moreover, due to the neo-liberal offensive (1980-1990) to the resumption of the 
"Fhcory of Human Capital", backed by Gary Becker, the wage differences came to 
be seen as responsibility of the worker himself, who thus became the chief culprit 
for all economic ills. The spreading of innovations due to technological changes (the 
convergence of information technology and telecommunications and the continuing 
cost reductions) greatly affected the rate of hiring as a function of the resizing of 
the human effectives face to an accelerated automation. Also, Security in all the 
countries that care about their retirees do not support the significant increase in 
life expectancy, discouraging retirements and thus retaining jobs occupied for much 
longer. All these factors contributed to a more chaotic work-employment scenario 
in Brazil from mid-1990 until about 2002 and they corroborate the interpretation of 
the plots. Many simulations are performed with vector V describing the geometry 
of the dynamics (theoretical) , intrinsic to the nature of the system itself, and vector 
U storaging the measurable content (observational) of the dynamics. For the more 
realistic Holling type II functional response, with fc 2 = 0, figures 2, 3, 4 and 5 show 
graphics of the general propensity of the system to less chaotic states accordingly 
the values 1.0, 1.4, 1.8 and 2.0 of the population interaction parameter a 2 - Figure 
5-a shows similar patterns for k 2 — 0.00008. Note that in all these figures there is 
a transition state to the phase of greater relaxation of the system. The great phase 
transition occurs close to the middle of the time series, from which the system goes 
on a relaxation stage. Comparing figures 1 and 5, the best fit realy seems to be 
for a 2 — 2.0 with phase transition around the zone between 50 and 100 months in 
figure 1. Figures 6, 7, and 8 show the system's phase space for a 2 = 1-0. The exact 
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fitting of this parameter is a matter of more studies, since predator-prey stability 
could come from many factors at play, such as migration of manpower, employers 
self-limiting their staffs, seasonal recruitment due to foreseen peaks of sales, and so 
on, in the same way as pointed out by Andrzej Pekalski^ln biological systems. 

5. Conclusion 

This study identified a predator-prey pattern in the work-employment system. The 
central equations applied are modified versions of those studied by Blasius et al, 
introducing pairwise dissimilarities between observations as a way to identify cyclic 
behaviors in the system. For carrying capacity ki = 0, simulations showed trajecto- 
ries with initial chaotic behavior and tendency to final states of greater equilibrium, 
which means that, even with zero carrying capacity of employers, the dynamics of 
the work-employment system in Brazil has demonstrated, during the past 12-15 
years, propensity to adjusting to less chaotic regimes. In part, this is due to the 
great migration in recent years to solidary economy and autonomous work. Present 
article analyzed the system dynamics under the interpretation of the population 
interaction parameter a,i as an index of the convergence between the aims of em- 
ployers and workers or as an index of the equilibrium created by lower demand for 
formal employment. Of course, interactions and proliferation of workers and em- 
ployers are not easily controlled, but the captivation of manpower is a variable that 
may be monitored and partially governed by public policies. Moreover, as Grafton 
and Echenique ^ pointed out, we have the problem to choose the more appropriate 
ecological model to describe the situation in study. In our opinion, as more is learned 
about the system, we must change adaptively the model in order to lower uncer- 
tainty. Far from exhausting the subject, it is clear that in its general form the model 
represents reasonably well the dynamics of predator-prey for the work-employment 
system in Brazil. We hope this work will serve as inspiration for further studies in 
order to compare the application of the model in different countries. 
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